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Abstract. Using the Becke-3-LYP functional, we have performed band structure calcu- 
lations on the high temperature superconductor parent compound, La2Cu04. Under the 
restricted spin formalism (p| = the R-B3LYP band structure agrees well with the 
standard LDA band structure. It is metallic with a single Cu — y^/O Pcr band crossing 
the Fermi level. Under the unrestricted spin formalism [p^ ^ p^), the U-B3LYP band 
structure has a spin polarized antiferromagnetic solution with a band gap of 2.0 eV, agree- 
ing well with experiment. This state is 1.0 eV (per formula unit) /otwer than that calculated 
from the R-B3LYP. The apparent high energy of the spin restricted state is attributed to 
an overestimate of on-site Coulomb repulsion which is corrected in the unrestricted spin 
calculations. The stabilization of the total energy with spin polarization arises primarily 
from the stabilization of the — y^ band, such that the character of the eigenstates at 
the top of the valence band in the antiferromagnetic state becomes a strong mixture of 
Cu — y^ jO Pa and Cu jO' pz. Since the Hohenberg-Kohn theorem requires the spin 
restricted and spin unrestricted calculations to give identical ground state energies and 
total spatial densities for the exact functionals, this large disparity in energy reflects the 
inadequacy of current functionals for describing the cuprates. This calls into question the 
use of band structures based on current restricted spin density functionals (including LDA) 
as a basis for single band theories of superconductivity in these materials. 



INTRODUCTION 



Almost immediately following the discovery of the superconducting cuprates (e.g. 
Lai_a;Baa;Cu04 in 1986), several research groups characterized the band structures within 
the density functional (DFT) formalism. ■'^"^ Using the standard local density approxi- 
mation (LDA), the resulting band structures consistently showed the Fermi level behav- 
ior of these materials to be characterized by a single metallic 2-D band, comprised of 
Cu — /O Pa hybrid orbitals (from here on called the — band). All other oc- 
cupied bands were buried 0.5 eV or more below the Fermi level. While the basic orbital 
picture appeared to agree with some experimental data, the absence of antiferromagnetic 
(AF) order in the band structure of the undoped parent compound i^e.g. La2Cu04) was 
cause for concern. In time, the inability of this band structure to explain an increasingly 
diverse range of normal state phenomena of the doped materials led many to conclude 
that Fermi liquid theory is not applicable. Indeed, no major theories based solely on 
this conventional band structure have survived. Still, the LDA band structure is widely 
viewed as a reasonable starting point for superconductivity theories and experimental data 
is routinely compared to these calculations. 

In the early 1990's, a number of groups succeeded in computing an alternative AF 
band structure for undoped La2Cu04 within the DFT formalism. Using a pseudopotential 
approach, Shiraishi, et al."^ achieved a spin polarized solution where j and | spins reside 
on different Cu sites of a doubled unit cell. This opened up a band gap of 0.6 eV between 
occupied and unoccupied bands (the measured gap is 2.0 eV^). Other groups attempted 
to correct certain known flaws in the LDA functional. Most notably, Svane^ applied a 
self-interaction correction local spin density (SIC-LSD) approach, in which the residual 
Coulomb interaction an electron improperly sees with itself is removed from the LSD 
functional. Spin localization was achieved with this method and an AF band structure 
was found with an indirect band gap of 1.04 eV. Temmerman, et al7 later found similar 
results using an alternative SIC-LSD approach. They reported an improved band gap of 2.1 
eV. Czyzyk and Sawatzky^ took yet another approach (LSDA-I-U), embedding a Hubbard 

2 



Hamiltonian into the Kohn-Sham LDA equation. They also achieved an AF state with a 
band gap of 1.65 eV. 

A common characteristic of all these calculations was a significant change in orbital 
character near the top of the valence band as compared to the standard LDA band struc- 
ture. All of the above authors noted a large increase in either the apical oxygen {O') or 
Cu density of states. While these results suggested the single band (x^ — j/^) picture of 
LDA may not be an adequate starting point for the doping range of superconductivity, it 
was not immediately apparent that a more complicated band picture was consistent with 
the experimental data either. Unresolved was the difficult question of how to describe the 
doped state of the superconductor which appears to produce a Fermi surface in the Bril- 
louin zone of the single unit cell.^ Removing electrons from a rigid band structure may be 
appropriate with the standard metallic state, but this procedure is less clear when starting 
from the undoped spin polarized insulating band structure in the reduced Brillouin zone 
(doubled unit cell). Thus, the LDA band structure has remained the de facto standard in 
the field to this day. 

In this work, we revisit the DFT band structure using the Becke-3-Lee- Yang-Parr 
(B3LYP) functional. The superiority of this hybrid functional, which includes a con- 
tribution of the exact Hartree-Fock exchange, has been well documented for molecular 
systems. Notably, Martin and lUas^^ showed the utility of such hybrid functionals in 
dramatically improving the calculation of the coupling constant J from cluster models 
of La2Cu04. Unfortunately, the method is still little used to determine band structures, 
partly because of the expense of such computations. We show here that spin unrestricted 
B3LYP (U-B3LYP) leads to an AF band structure in agreement with the SIC-LSD and 
LSDA+U results cited above. This assuages doubts about that work and confirms the 
Cu jO' Pz character at the Fermi level. Most importantly, this AF state is found to 
be 1.0 eV per formula unit more stable than the state calculated from the spin restricted 
(R-B3LYP) functional. Considering the Hohenberg-Kohn theorem^^ requires the exact 
spin restricted and spin unrestricted functionals to yield identical total spatial densities 
{p^ + Pl) and total energies, this discrepancy indicates a serious flaw in current functionals 
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(B3LYP and LDA). Therefore, use of the LDA band structure to justify single band models 
of superconductivity in the cuprates is highly questionable. 

RESULTS AND DISCUSSION 

Calculations were performed using CRYSTAL98^^ which employs an atomic Gaussian 
type orbital basis set. For O, the standard 8-411G basis set with a D polarization exponent 
of 0.65 was used.^^ For Cu and La, the Hay and Wadt^^ effective core potentials (ECP's) 
were used. These ECP's treat explicitly the outer core (3s and 3p for Cu, 5s and 5p for 
La) and valence electrons. The basis sets used with these ECP's were modified from the 
original basis sets of Hay and Wadt, since some functions are too diffuse for calculations on 
crystals. For Cu, the two diffuse S exponents were replaced by a single exponent optimized 
to 0.30 from LDA calculations on La2Cu04. The two Cu diffuse P exponents were replaced 
by a single exponent optimized to 0.20. The basis set was contracted to (3s3p3d) based on 
atomic Cu(n) calculations. For La, the two diffuse S exponents were replaced with a single 
exponent optimized to 0.10. The two diffuse P exponents and the diffuse D exponent were 
removed without replacement. The basis set was contracted to (3s2pld) based on atomic 
La(ni) and La(II) calculations. Overall, the quality of the basis set is superior to that 
used by Su, et al}^ in their CRYSTAL95 Hartree-Fock (HF) study of La2Cu04. Several 
alternative basis sets were tested, all leading to similar results. The tetragonal La2Cu04 
crystal structure was taken from Hazen.^''' 

Figure 1 presents the results of our LDA and R-B3LYP calculations with restricted 
spin and tetragonal symmetry. The LDA band structure is in excellent agreement with 
previous plane wave calculations^"^ and there is little difference with the R-B3LYP band 
structure result. As expected for restricted spin calculations, both methods produce Pauli 
paramagnetic (PM) band structures. In each, the only band crossing the Fermi level is the 
highly 2-D — y^ band. The next band (^^) is approximately 1 eV below the Fermi level. 

The lack of antiferromagnetic order in these band structures limits their usefulness for 
analyzing properties of the undoped material. On the other hand, such PM band structures 
may be entirely appropriate for understanding the nature of the doped state. Indeed, 
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the LDA and R-B3LYP band structures seem to be consistent with some experimental 
data such as the ARPES Fermi surface,^ but they remain incompatible with many other 
experiments. An analysis of the density of states shows the nature of the intrinsic undoped 
hole (lowest unoccupied states totaling 1 hole per formula unit) is 48% Cu — y"^ and 
47% O Pfj. While we have not explicitly carried out computations on the doped state 
(a;=0.15), we can estimate the nature of the doped hole using a rigid band model. The 
character of this doped hole (highest occupied states totaling 0.15 electron) is 47% Cu — 
y^, 38% O Pa, 6% O' pz, and 5% Cu z^. The relatively small amount of O' and 
character in these orbitals is in keeping with most models for superconductivity in the 
cuprates. However, x-ray absorption studies (XAS) support a total z'^ hole contribution 
of 5% to 20% and a similar range for O' Pz-^^ Details such as the ARPES pseudogap/^ 
the anomalous background signal, ^° and other probes of the normal state properties such 
as the NMR,^-*^ resistivity,^^ and neutron scattering^^ also appear to have no explanation 
using this conventional band structure. 

We find that, a diff'erent state emerges under the unrestricted spin (U-B3LYP) imple- 
mentation of this functional. Using an orthorhombic unit cell, the resulting band structure 
from these calculations is shown in Figure 2. The combination of a spin functional and 
a doubled unit cell allows for a possible spin polarized solution. Indeed, we find an ap- 
parent AF state with a gap of 2.0 eV. The band dispersion is in excellent agreement with 
previously published DFT band structures for this AF state, and the computed gap 
agrees with the measured gap.^ Furthermore, a U-B3LYP calculation of the ferromag- 
netic (FM) state, which has a pure t spin per formula unit under tetragonal symmetry, is 
found to be 0.18 eV higher in energy. This compares favorably (noting caveats) with the 
experimental J value of 0.13 eV.^^ 

A principal advantage of this new U-B3LYP AF band structure is that it follows un- 
ambiguously using a well established functional. No additional empirical corrections were 
necessary. Furthermore, while no comparison of the relative stabilities of the LDA state 
and the SIC-LSD (or LSDA+U) state has been previously reported, such a comparison 
is rather straightforward with our calculations. Significantly, the U-B3LYP AF ground 
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state and FM excited state are found to be 1.0 and 0.82 eV per La2Cu04 formula unit 
more stable than the R-B3LYP state. This represents a rather dramatic failure on the 
part of the spin restricted functional. Clearly the U-B3LYP calculation leads to a superior 
representation of the ground state of La2Cu04. 

To better understand the significance of this 1.0 eV energy difference, one must 
consider the theoretical foundation of DFT. The Hohenberg-Kohn theorem^° proves the 
existence of two functionals Fhk[p] and FHK[pt'Pi] where the first is a functional of the 
total density p = Pt + Pj. and the second is a functional of the two spin densities p-|- and p|. 
For a spin-independent potential, the first functional will lead to the exact ground state 
energy and total spatial density, p = p-f+pj,. For spin-dependent spatial potentials, the use 
of the second functional will lead to the ground state energy and spin densities p-|- and pj^. 
In addition, FHK[pt)Pi] can be used for the special case of a spin- independent potential 
where it must obtain the same ground state energy and total density p as Fhk[p]- 

The Kohn-Sham orbitals arising from energy minimization for La2Cu04 will be of 
restricted (doubly occupied) Hartree-Fock type for Fhk[p] with p(r) = 2 X^occ 
of unrestricted Hartree-Fock type for FhkIpT) Pj.] where the t spin orbitals may be different 
from the J, spin orbitals. Both minimizations should lead to exactly the same energy and 
ground state total density p. If the R-B3LYP functional were close to Fhk[p] and the 
U-B3LYP functional were close to Fhk[pt;Pj.]; then the energies of the two calculations 
for La2Cu04 should also be close. 

Our computed 1.0 eV energy difference leads us to conclude that at least one of the 
two functionals is not close to the exact functional. Since the U-B3LYP band structure is 
an excellent description of the AF state, we are led to question the quality of the R-B3LYP 
functional. 

It is well known in fact that spin density functional approximations to Fhk [pt ; Pi\ are 
superior to total density approximations to Fhk[p] since it is much easier to empirically 
design functionals to correct for the exchange coupling of like spins when p is separated into 
p-|- and p|. In the restricted spin functionals, R-B3LYP and LDA, the improper treatment 
of exchange in strongly correlated systems leads to an overestimate of on-site Coulomb 
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repulsion. 

A classic example of the problems that occur with these functionals is the dissociation 
of H2. At equilibrium distances, the molecule is well described both by closed-shell LDA 
and by B3LYP functionals. However, at the dissociation limit both functionals lead to an 
energy calculated to be higher than two H atoms due to the ionic components (H"*" + H~). 
With current functionals, to properly describe dissociation requires calculation of either 
the triplet state or the symmetry broken unrestricted "singlet" spin state, either of which 
leads to net formation of magnetic moments on each H atom. This difference between the 
spin restricted and unrestricted functionals should not be considered a failure of density 
functional theory. Instead, it is a failure of the empirical spin restricted functionals that 
have been developed so far. 

By analogy to H2, the FM state of undoped La2Cu04 (where each site has a pure t 
spin) and the AF state (where each site has either a pure t or a pure J, spin) are well de- 
scribed by the spin unrestricted U-B3LYP functional. Using the spin restricted formalism 
where each site is 50% j AND 50% | (R-B3LYP), leads not only to an overestimate (1.0 eV 
per formula unit) of the total energy, but also to an incorrect band structure. In particular, 
the x'^ — band is elevated relative to the other bands as a result of the improper on-site 
Coulomb repulsion associated with the Cu — y^ jO p^j orbital. Removing this repulsion 
through the localization of spins in either the FM or AF states stabilizes the x"^ —y^ band, 
resulting in an increase in character near the top of the valence band. 

As compared to the R-B3LYP state, the added stability of the U-B3LYP total energy 
by 1.0 eV is associated with an approximately 1 eV stabilization of the — band with 
respect to the other bands. This is most noticeable in the nature of the doped hole. We 
characterize the U-B3LYP undoped hole as 56% Cu x"^ — y^ and 38% O a picture 
not substantially different from the R-B3LYP calculation. The ratio of Cu a;^ — y^ /O Pa- 
character is somewhat larger in the U-B3LYP state, but qualitatively both calculations 
agree that the undoped hole states are nearly purely derived from these two orbitals. The 
picture changes substantially upon considering the nature of the doped holes. Figure 3 
presents a detail of the density of states in the vicinity of the Fermi level for the a;=0.15 
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doped state, assuming a rigid band model. The nature of the doped hole is characterized 
as 17% Cu - y^, 40% O p^, 21% O'vz, and 19% Cu . This is summarized in 
Table I. The significant increase in the Cu and O' pz character of the doped hole 
as compared to the R-B3LYP results is comparable to that noted in previous AF band 
structure calculations. While doubts about the quality of the previous band structures 
undermined the significance of these findings, the cumulative weight of these results now 
strongly favors the scenario where holes are formed upon doping. 

While the U-B3LYP band structure may be a good representation of the AF ground 
state of undoped La2Cu04, we should not lose sight of the fact that ultimately a spin 
restricted (total density) functional should be equally successful. Indeed, a proper PM 
band structure is more useful in understanding the role of doping in superconductivity. 
In our view, such a band structure has yet to be achieved with DFT, leading many to 
conclude that Fermi liquid theory has failed for these materials. The simpler answer might 
just be that LDA has failed. The calculations presented here demonstrate the shortcomings 
of these methods and suggest possible solutions. We might postulate the existence of a 
modified R-B3LYP functional which gives exactly the same total density and total energy 
as the U-B3LYP functional. The PM band structure associated with this new functional 
may be substantially different from the one we've calculated here. Elsewhere, we have 
already incorporated such a correction into a simple tight binding model for the doped 
superconductor.^^ By effectively introducing a local magnetic moment at each Cu site 
within a PM model of optimally doped Lai.85Sro.i5Cu04, we showed that the x'^ — 
band is significantly stabilized relative to the other bands. This brings the narrow ^ band 
to the Fermi level. The resulting band structure has a unique crossing between the 2-D-like 
— band and the 1-D-like z"^ band. This reflects a dramatic first-order correction to 
the standard band structure. Indeed, this model has already been used to interpret the 
ARPES pseudogap and anomalous background, the NMR Cu and O relaxations and 
Knight shifts, the Hall effect, and Josephson tunneling^® have already been presented 
with this model. 
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CONCLUSIONS 

In summary, we have presented the results of R-B3LYP and U-B3LYP band structure 
calculations on La2Cu04. The R-B3LYP results are in good agreement with previous 
LDA calculations and the U-B3LYP results are in good agreement with previous SIC- 
LSD calculations (among others). The large discrepancy in energy between the two states 
(1.0 eV per formula unit) is attributed to an improper overestimate of on-site Coulomb 
repulsion within the spin restricted calculations. The automatic correction of this error 
within the U-B3LYP AF state leads to the stabiliation of the x'^ — band relative to the 
other occupied bands. As a result, the band is then brought to the top of the valence 
band and contributes significantly to the doped hole states. These results cast significant 
doubt on the continued use of LDA band structures as the starting point for theories of 
superconductivity in these materials. 
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gene Heifets. This work was partially supported by the Materials and Process Simulation 
Center (MSC) at Caltech which is supported by grants from DOE- ASCI, ARO/DURIP, 
ARO/MURI, 3M, Beckman Institute, Seiko- Epson, Dow, Avery-Dennison, Kellogg, and 
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Table I. Orbital character of intrinsic undoped holes (totaling 1 hole), doped holes (total- 
ing 0.15 hole), and total holes at optimal doping (totaling 1.15 holes). Results are shown 
for both the restricted spin R-B3LYP state and unrestricted spin U-B3LYP AF state. 



Orbital 


undoped 


B3LYP 

doped 


total 


undoped 


UB3LYP 

doped 


total 


Cu - 


48% 


47% 


48% 


56% 


17% 


51% 




47% 


38% 


46% 


38% 


40% 


38% 


O' p. 


1% 


6% 


2% 


1% 


21% 


4% 


Cu 


1% 


5% 


2% 


0% 


19% 


3% 
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FIGURE CAPTIONS 

Figure 1. Band dispersions of La2Cu04 plotted along symmetry lines of the tetragonal 
Brillouin zone (see reference 1) from restricted spin (a) LDA and (b) R-B3LYP calculations. 
Results are in good agreement with the LDA computations of references 1-3. 
Figure 2. Band dispersion plotted along symmetry lines of the orthorhombic Brillouin 
zone (see reference 8) from the unrestricted spin U-B3LYP calculation of La2Cu04. Results 
are in good agreement with the SIC-LSD and LSDA+U computations of references 6-8. 
Figure 3. Detail of the density of states from the unrestricted spin U-B3LYP AF calcu- 
lation of La2-a;Sra;Cu04 (x = 0). Fermi level is positioned to a doping level of x=0.15. (a) 
projected Cu DOS vs. total DOS. (b) projected O DOS vs. total DOS. (c) projected O' 
DOS vs. total DOS. (d) projected Cu x'^ - y'^ DOS vs. total Cu DOS. (e) projected Cu 
z'^ DOS vs. total Cu DOS. 
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